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We compute cosmic ray (CR) nuclei, proton, antiproton, electron and positron spectra below 1 
TeV at Earth by means of a detailed transport description in the galaxy and in the solar system. CR 
spectra below 10 GeV are strongly modified by charge-sign dependent propagation effects. These 
depend on the polarity of the solar magnetic field and therefore vary with the solar cycle. The 
puzzling discrepancy between the low-energy positron fraction measured by PAMELA and AMS-01 
is then easily explained by their different data-taking epochs. We reproduce the observed spectra 
of CR light nuclei within the same galactic and solar-system propagation model. 
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Introduction: The propagation of Cosmic Rays (CRs) 
in the Galaxy is far from being fully understood. One 
major challenge is represented by the spectrum of the 
positron fraction (PF), e + /(e~ + e + ), measured now by 
PAMELA and Fermi with high accuracy at energies rang- 
ing from below 1 GeV up to about 100 GeV [IM], and 
previously by AMS-01 [5] and many other experiments. 
While the well known increase of the PF with increasing 
energy above 10 GeV has attracted most of the attention 
of the astroparticle physics community, the disagreement 
between PAMELA and AMS-01 PF below 10 GeV has 
motivated several investigations 6 9 but has never been 
convincingly explained so far. 

Before they are detected at Earth, CR e~ and e + lose 
energy due to solar winds while diffusing in the solar 
system [10]. This modulation effect depends, via drifts 
in the large scale gradients of the solar magnetic field 
(SMF), on the particle's charge including its sign [TT] . 
Therefore, it depends on the polarity of the SMF, which 
changes periodically every ~11 years [12] . Besides the 
11 year reversals, the SMF has also opposite polarities in 
the northern and southern hemispheres: at the interface 
between opposite polarity regions, where the intensity of 
the SMF is null, a heliospheric current sheet (HCS) is 
formed (see e.g. [12] )• The HCS swings then in a region 
whose angular extension is described phenomenologically 
by the tilt angle a. The magnitude of a depends on solar 
activity. Since particles crossing the HCS suffer from 
additional drifts because of the different orientation of 
the magnetic field lines, the intensity of the modulation 
depends on the extension of the HCS. The PF spectra 
measured by AMS-01 and PAMELA can therefore differ 
because they were measured in different periods of the 
solar activity [TllSlfTl]. 

We show in Fig. [T] how changing polarity affects the 
PF, computed using the methods that will be described 
later. These effects are clearly evident below 10 GeV in 
the model curves and are large enough to explain the 
discrepancy between AMS-01 and PAMELA data, which 
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FIG. 1. The positron fraction measured by Fermi, PAMELA 
and AMS-01 is shown. The LIS is shown as the red dashed 
curve. Solid curves show the Earth positron fraction com- 
puted evolving the LIS for a — 30° and positive polarity 
(black) or negative polarity (violet). The long-dashed blue 
curve represents a prediction for AMS-02, which is taking 
data in a period of negative polarity and high solar activity. 



were taken in periods of opposite polarity and compara- 
ble tilt angle a ~ 30°. Indeed, we reproduce very well the 
PF measured by AMS-01 with the solar model adapted 
to AMS-01 data taking conditions, but remarkably we 
also nicely reproduce PAMELA data by tuning the solar 
model to the appropriate PAMELA period. This con- 
firms that the low energy PF can be interpreted self- 
consistcntly by taking into account drift effects in solar 
propagation. We also show the PF computed for a = 75° , 
appropriate for year 2012, during which AMS-02 has been 
taking data [T5]. This prediction will soon be tested by 
that experiment. 

This result is the main achievement of this Letter. We 
will support it by describing the computational methods 
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and assumptions we adopted to compute the PF. The 
main novelty with respect to previous investigations [HE] 
is that for the first time we have computed the PF within 
a galactic+solar propagation model self-consistently de- 
rived from observations of CR nuclear and lepton spectra. 
Our results demonstrate that having a precise account of 
solar effects, going beyond the force-field model of [TU] . 
is crucial to interpret the data available nowadays. 

Method: We first select one propagation model by im- 
posing that CR nuclear spectra at Earth be correctly 
reproduced. We compute the galactic interstellar spec- 
tra (LIS) first and then apply a detailed 4D (3 spatial + 
energy) solar propagation model to CR nuclei to obtain 
the Earth spectrum to be compared with observations. 
Then, within the propagation model determined in this 
way, we compute the e + and e~ spectra and the PF. 

Both propagation regimes are described by a gen- 
eral diffusion-convection-reacceleration-energy-loss equa- 
tion [16l [17]. We propagate CR species in the Galaxy 
starting from Z = 28 with the DRAGON code [18]. To 
avoid boundary effects, we set the vertical boundary of 
the numerical domain L = 2zt, with z t the half-height 
of the diffusion region. We assume a diffusion coefficient 
in the Galaxy D(r, z,p) = D /3 v (p/po) S exp(z/z t ), where 
p(p) is particle rigidity, /3 is particle's velocity, po = 3 GV, 
77 and S are constant parameters. For reacceleration, we 
use the standard description in terms of momentum dif- 
fusion, with D pp cx p 2 v 2 A /D, where va is the Alfven ve- 
locity. Convective effects due to stellar winds are de- 
scribed by a convection velocity directed along z outside 
the galactic plane vc(z) — vco + wciG^I/l kpc). 

We set S = 0.5, z t = 4 kpc, D = 2.5 x 10 28 cm 2 /s, 
va = 18 km/s, vco = and vci = 5 km/s/kpc, 77 = 1. 
We assume for all nuclei a power-law injection spectrum 
with spectral index 7 = 2.3. For primary electrons, 
we assume a broken power-law with 7 = 1.5/2.65 be- 
low/above 6 GeV. The presence of such a break in the 
primary e~ spectrum is also confirmed by the analysis of 
the diffuse galactic synchrotron emission [H)H2"l"] . In or- 
der to reproduce high energy data, we also consider an 
extra-component of primary electrons and positrons with 
7cxtra = 1.5 and exponential cutoff at 1.2 TeV [22] 123], 
We fix its normalization so that we match the Fermi 
e + + er spectrum at 300 GeV, while the normalization 
of the primary electron spectrum is fixed at 33 GeV. 

For the 4D propagation of CRs in the solar system we 
develop instead a new numeric program, HelioProp. 
We follow the stochastic differential equation approach 
of [2~4Tf2"6"] (see [23 [25] for a general description) , where 
the CR phase-space density is computed by sampling and 
averaging upon pseudo-particle trajectories. Each tra- 
jectory is in fact the result of a deterministic component 
related to the drifts, and of a random walk component, 
whose amplitude is sampled according to the local diffu- 
sion tensor. Pseudo-particles injected at the Earth posi- 
tion are followed backward in time during their propaga- 



tion in the solar system until they reach the heliopause, 
where their properties are recorded (we refer the reader to 
[24TI26] for more details on the actual numerical scheme 
to be used). We propagate 10 4 pseudo-particles per en- 
ergy bin, which allows us to efficiently sample the Green 
function and to keep statistical fluctuations below 1% 
[2"4l |2"5"] . The LIS flux, which is effectively a boundary 
condition for this problem, is then used as an appropri- 
ate weight to determine the Earth spectrum. 

We specify our model for solar propagation by fix- 
ing the solar system geometry, the properties of diffu- 
sion and those of winds and drifts. For simplicity we 
assume spherical geometry, although some latitudinal de- 
pendence of the position of the heliopause, which we set 
at 100 AU, has been experimentally found (see [29] and 
Refs. therein). 

We describe the solar system diffusion tensor by 
K(p) = diag(if || , K± r , K±g)(p), where || and _L are set 
with respect to the direction of the local magnetic field. 
We assume no diffusion in the _L tp direction and we de- 
scribe as drifts the effect of possible antisymmetric com- 
ponents in K. For the parallel CR mean- free-path we 
take A,, = A (p/1 GeY)(B/B(g)~\ with A = 0.15 AU 
and £?0 = 5 nT the value of the magnetic field at Earth 
position, according to [24l [25]. For p < 0.1 GeV, X\\ 
does not depend on rigidity. The value of Ao and the 
rigidity dependence of Ay are compatible both with the 
measured e~ mean- free-path (see, e.g., [3U]) and with 
the proton mean-free-path inferred from neutron moni- 
tor counts and the solar spot number [29] . We then com- 
pute if II = \nv/3. Perpendicular diffusion is assumed 
to be isotropic. According to numerical simulations, we 
assume \± r ,e = 0.02A|| [3"T] . 

For the SMF, we assume a Parker spiral, although more 
complex geometries might be more appropriate for peri- 
ods of intense activity 

f r V 2 A ftrsintf A 

where fl is the solar differential rotation rate, 8 is the 
colatitude, Bq is a normalization constant such that 
\B\(1 AU) = 5 nT and A = ±H(0 - 9') determines 
the MF polarity through the ± sign. The presence of 
a HCS is taken into account in the Heaviside function 
H(9 — 6'). The HCS angular extent is described by the 
function 8' = 7r/2 + sin~ 1 (sinasin(<y9 + Qr/Vsw)), where 
< a < 90° is the tilt angle. The drift processes occur- 
ring due to magnetic irregularities and to the HCS are re- 
lated to the antisymmetric part Ka of the diffusion tensor 
as m "drift = V x (K A B/\B\) = sign(<7>/3V x (r L s), 
where Ka — pv /iqB 1 tl is the particle's Larmor radius 
and q is its charge. We refer to [24l|25] for more details on 
the implementation of the HCS and of drifts. Adiabatic 
energy losses due to the solar wind expanding radially at 
Vsw ~ 400 km/s are taken into account. 
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FIG. 2. Time evolution of the tilt angle. Vertical bands high- 
light data taking periods of relevant experiments. Data rele- 
vant for the positron fraction were taken by PAMELA in the 
period highlighted by the red area. The larger greyish shad- 
owed area marks the data-taking period relevant for PAMELA 
absolute e~ spectrum data. 



We report in Fig. [2] the values of a [12] inferred from 
solar models, together with the periods in which rele- 
vant CR experiments were taking data. Predictions from 
the "L" (line-of-sight boundary conditions) model seem 
more accurate for periods of decreasing tilt angle, while 
for periods of increasing a model "R" (radial boundary 
conditions) is more precise (33J [34] . 

A and a are of particular importance for CR propaga- 
tion in the heliosphere. If q-A < 0, drifts force CRs to dif- 
fuse in the region close to the HCS, which enhances their 
effective propagation time and therefore energy losses, 
while if q ■ A > drifts pull CRs outside the HCS, where 
they can diffuse faster [211 [22] • 

For each CR experiment we fix the relevant solar prop- 
agation model as follows: a = 30° for both PAMELA 
and AMS-01 data, but A = + for AMS-01 and A = - 
for PAMELA. For ACE [35] we use a = 10° and A = +. 
We do not vary the intensity of the magnetic field at 
Earth position. 

Results: We obtain a good description of the B/C spec- 
trum over the entire energy range, also including ACE 
data [35] , as we show in Fig. [3] We achieve this without 
resorting to large va, which is in conflict with the diffuse 
synchrotron galactic emission [201 121] > nor to modified 
low-energy diffusion. Combinations of these possibilities 
have been invoked to reproduce ACE data [HrE]. We 
checked that also the spectra of C, O and 10 Be/ 9 Be are 
correctly reproduced. 

The proton and He spectra measured by PAMELA are 
reproduced very well, without introducing a break in the 
injection spectrum at ~ 10 GeV and with the same injec- 
tion index, as shown in Fig. [4] Interestingly, the AMS-01 
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FIG. 3. Experimental data for the B/C ratio [36H39] com- 
pared to our model. 
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FIG. 4. Experimental data for proton spectra from PAMELA 
and AMS-01 and for the He spectrum from PAMELA com- 
pared to theoretical calculations for a = 30° . 



and PAMELA proton spectra are very similar at low en- 
ergy, although they have been taken in different periods 
of solar activity. Given the observed values of a for the 
respective data-taking periods, we cannot reproduce both 
measurements with our model. To achieve a good match 
of both, we would need to assume a < 5° for AMS-01, 
which is not inconsistent with observations given that the 
AMS-01 data taking period followed a phase of very low 
solar activity and that the HCS takes about 14 months 
to propagate from the Sun to the Heliopause 29 . This 
has however negligible impact on the computation of sec- 
ondary e + and e~ in the energy range relevant for the PF, 
because the discrepancies are at E < 3 GeV. BESS spec- 
tra |40| . which have been taken over a period spanning 
almost a complete solar half cycle, might help under- 
stand this issue, but given that we are interested here in 
the comparison of the PF of PAMELA and AMS-01, we 
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FIG. 5. Experimental data for the p/p ratio from PAMELA 
and BESS [40] compared to theoretical calculations. In order 
to highlight charge-sign effects, models for opposite magnetic 
polarity and same a = 30° are shown. 
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FIG. 6. The absolute e and e spectra measured by different 
experiments are compared with our calculations. 



prefer to leave the study of BESS spectra for future work 
and compare only PAMELA and AMS-01 data. Another 
possibility is to more finely tune the propagation model. 
Indeed, by slightly adapting the propagation parameters 
a much better fit of both PAMELA and AMS-01 proton 
data can be achieved, at the price of a slight worsening 
of the B/C. However, we explicitly verified that our fit of 
the PF is not significantly affected. 

Also p (not shown) and p/p spectra are very well repro- 
duced (see Fig. [5| . In this case the separation between 
model calculations in different periods of solar activity is 
comparable to experimental errors. 

We show finally in Fig. [6] the absolute e~ and e + spec- 
tra. We achieve a rather good agreement, within experi- 
mental uncertainties, with the observed absolute e + spec- 



trum, which is however slightly underestimated between 
1 and 3 GeV. The e~ spectrum below 5 GeV is under- 
estimated. In order to better understand this issue, we 
consider a modified model, in which we lower z t = 2 kpc 
and set 77 = —0.2. In order to ensure that the B/C ra- 
tio is still reproduced, we keep D /z t constant [IS] and 
slightly readjust va — 13 km/s. With these settings, we 
still correctly reproduce the B/C, p and p data and affect 
only marginally the PF, but we improve the e + and the 
e~ spectra. This may hint at low-energy modified diffu- 
sion effects being indeed relevant, but we must remark 
that theoretical uncertainties on secondary e fluxes are 
large in this energy range because of poor knowledge of 
the production cross sections in pp collisions [6] [9] . 

Conclusions: Accounting for CR propagation in the 
solar system beyond the force-field description, is very 
important to understand present-day observations. We 
have computed the local interstellar spectra of several 
CR species and propagated them in the heliosphere by 
accounting for charge-sign dependent drifts. We applied 
our solar model to the B /C for the first time and showed 
that ACE observations can be interpreted within galactic 
propagation models in which also the observed spectra 
of p, He and p can be successfully reproduced without 
requiring spectral breaks. 

We then computed the CR e~ and e + spectra and the 
positron fraction within the same models. The discor- 
dant observations at low energy by PAMELA and AMS- 
01 can be explained by charge-sign dependent drifts. For 
the first time, this result is achieved within galactic prop- 
agation models which are self-consistently derived from 
observations of CR nuclei, and by solely adjusting the 
solar propagation model to match the actual observation 
conditions of each experiment. 

The e~ and e + low energy spectra are not satisfacto- 
rily reproduced, unless modified low-energy diffusion and 
a rather small halo are assumed. Enhancing the produc- 
tion of secondary leptons might improve the e + spectra, 
but has an impact on the e~ spectrum only at very low 
energies. This might hint at the need of a more elabo- 
rated description of galactic and solar propagation, or of 
a more detailed analysis of secondary production cross 
sections. High quality and high statistics AMS-02 [TS] 
data, which are taken in a period of high solar activity 
and changing polarity, will help clarifying these issues. 

We have shown results only for one particular set of 
propagation parameters. In fact, we tested several com- 
binations of solar and galactic parameters, finding that 
in many cases a good agreement with all datasets can be 
achieved. A more detailed study of the parameter space 
is left for future work. 
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